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We are interested in the stability of single localized structures in a real Swift-Hohenberg equation 



subjected to a delayed feedback. We shall show that variation in the product of the delay time and 
jyy \ the feedback strength leads to complex dynamical behavior of the system, including formation of 

<jj ; Solitm ringS , ^ ^ OT _ We _ d e . ^ ^ 

of the delayed system and derive a system of order parameter equations for the position of the 
localized structure as well as for its shape. In a special case, a normal form of the delay-induced 
*i> ' drift-bifurcation is obtained, showing that spontaneous motion to the lowest order arises without 

V) ■ change of shape. 
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Control and engineering of complex spatio-temporal patterns in high- dimensional non- 

nn 

equilibrium systems has evolved as one of the central issues in applied nonlinear science |]J, |2|] . 
Among other control techniques, a quite simple and efficient scheme is the time delayed 
feedback (TDF) [3]. Although TDF method was originally designed to control dynamical 
systems with a few degrees of freedom, it has been successfully applied to a large number of 
theoretical and experimental high-dimensional spatially extended systems Q). In particular, 
control of dissipative localized structures (LSs) has been of increasing interest in recent 
years 0, Q|. Among other LSs, traveling pulses and fronts subjected to a TDF appear in 
various contexts [3] . Recently, properties of two-dimensional LSs by means of a real Swift- 
Hohenberg equation subjected to TDF were studied in [8[. It was shown that when the 
product of the delay time and the feedback rate exceeds some critical value, a single LS 
starts to move in an arbitrary direction. Moreover, an analytical formula for its velocity was 
derived. We refer to such kind of instability as delay-induced drift-bifurcation. 

In this Letter we investigate the stability of a single LS in a real Swift-Hohenberg equation 
subjected to a TDF in details. Varying of the product of the delay time and the feedback 
rate, among moving LSs, also leads to the formation of complex spatio-temporal patterns, 
such as oscillons, soliton rings or labyrinths. Moreover, we provide a bifurcation analysis of 
the delayed system and derive a set of order parameter equations for the position of the LS as 
well as for its shape. In a special case, a normal form of the delay-induced drift-bifurcation 
is also derived, showing that the drift-induced motion can arise without pronounced change 
of the LS's shape. The bifurcation analysis is performed in general form and therefore can 
be applied to any system possessing solutions in form of LSs. 



We start with the delayed Swift-Hohenberg equation (DSHE) js] 

— = -01 V 2 X - a 2 V 4 X + f(X) + a (X(t) - X(t - r)) , (1) 

where X := X(x, £), x = (x, y) T is the real distributed order parameter, describing the 
deviation of the cavity field, ai, a 2 > 0, r is the delay time, whereas a denotes the positive 
delay rate. In the polynomial nonlinear function f(X) = Y + C X — X s the coefficient 
Y describes the deviation of the injection field and C is the cooperative parameter. The 
derivation of the Eq. (JB as well as detailed description of physical meaning of the param- 
eters can be found in p, [gj]. For a = Eq. (pQ) becomes the Swift-Hohenberg equation 
(SHE), describing nascent optical bistability in the weak dispersion limit [10]. Note that the 



SHE possesses a Lyapunov functional, i.e., no Hopf bifurcations as well as time-dependent 
solutions are possible. Furthermore, stable stationary LSs are homoclinic solutions of SHE 
with dX/dt = 0. Notice that in the presence of TDF Eq. (pQ) loses the gradient structure. 
Although the homogeneous steady states as well as stationary localized solutions are not 
affected by the TDF, their stability properties may change. 

In order to investigate the stability of both stationary and homogeneous solutions we rewrite 
Eq. (pQ) in the general form for the vector function q = q(x, t),x= (x, y) T 

&q(x, t)=£ q(x, t) + a E (q(x, t) - q(x, t - r)) . (2) 

Here, £ denotes the nonlinear Swift-Hohenberg operator and E is the identity matrix. As- 
sume that the stationary solution of (Ej) qo(x) exists for all values of a. In the simplest case 
it is a single stationary LS with rotational symmetry. For a = linear stability analysis 
of q yields a real spectrum, given by the linear eigenvalue problem £'(qo) V = I^^P-> since 
the linearization operator £'(qo) is self- adjoint. Due to translational invariance \i = is 
an eigenvalue of £'(qo), corresponding to two independent neutral eigenfunctions cp^, also 
referred to as Goldstone modes. The discrete spectrum of £'(qo) in two dimensions contains 
modes oc <£> n (x) e m< ^ with n = 0, ±2, where the mode with n = results in the change of the 
size of the LS, and n = ±2 lead to deformations. All other modes are strongly damped. We 
suppose that for r = the stationary solution q is stable. For the spectrum of the 

corresponding linear operator can be constructed from the spectrum of £'(qo)- Indeed, both 
operators possess the same set of eigenfunctions v?( x )> whereas the complex eigenvalues A 
of the delayed problem are determined as solutions of the transcendental equation [8[: 

/x = A - a (1 - e~ Xr ) , (3) 

n i / / 

that can be found in terms of Lambert W-functions [11] as A = /Li+a+^ W m [—ar exp(— 
a)r)) , m G Z. Hence, the stationary solution of © remains stable, if Re(A(/x)) < for 
any \i from the spectrum of £'(q ). Bifurcation points correspond to /x, at which Re(A(/i)) 
vanishes. In particular, instability criterion for all eigenvalues A, corresponding to the real 
eigenvalues /i, is given by the relation a := ar > max{— fir/2 , l} . One can see that 
H = yields the bifurcation point a = 1, corresponding to the onset of spontaneous motion, 
observed in [8]. Furthermore, eigenvalues ji of discrete critical modes n = 0, n = ±2 are 
close to zero, i.e., for typical system parameters the instability threshold yields the value 
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Figure 1. A real part of the dispersion function A(fc), calculated at a = 1.05 for three different values 
of the delay time r, corresponding to a homogeneous Hopf bifurcation with k = for r = 0.02 
(dashed rot line), a traveling wave bifurcation with k = k c for r = 0.048 (dotted blue line) as well 
as stable solution qh for r = 0.24 (black dash-dotted line). Only the main branch of Re(A(fc)) is 
shown. A gray solid line corresponds to the function ji(k) for r = 0. 

a — 1, which becomes the bifurcation point of higher codimension. In addition, using a as a 
control parameter, one can find a critical delay time r c := xaj \i for which each of the modes 
becomes unstable by means of the solvability condition: 

± arccos(l + x) + 2nm = ±a y/l — (1 + x) 2 , m e Z . (4) 

The next point to emphasize is that the TDF can also affect the stability of the homogeneous 
solution qh of the system (J2D [9]. Suppose that qh is stable for r = 0, i.e., the dispersion 
function ji{k) = a\k 2 — a 2 k 4 + /'(qh) < for all values of k. For r > the linear stability 
analysis with respect to perturbations oc exp(i/cx + A(fc)t) yields the same Eq. ([3]). The 
dependence of Re(A(fc)) on the delay time r, calculated for fixed value of a = 1.05 is 
presented in Fig. [TJ One can see that for r = qh is stable. However, for r > the behavior 
of Re(A(fc)) indicates the existence of a pair of homogeneous Hopf bifurcations (HH) with 
k = (see Fig. Q] for r = 0.02) as well as traveling wave bifurcations (TW) with a finite 
wave number k = k c = ±^/ai/2a 2 (r = 0.048) [9]. The thresholds, associated with these 
instabilities can be found within the solvability condition (|4]) . An increasing of r leads to the 
stabilization of qh (r = 0.24). Notice that the behavior of Re(A(fc)) for large k is different 
for different values of r. While for small delays Re(A(fc)) decays as fc 4 , it decays rather as k 2 
for larger r, which is more in the nature of e.g., reaction-diffusion systems. Hence, in this 
parameter region one would expect solutions which are not generic for the SHE. 

In Fig. [2] a bifurcation diagram in (a, r) plane is presented. Each line on the plot is a 

4 



1 5 



10 15 20 

a 



Figure 2. Bifurcation diagram in (a, r) plane. Different lines separate instability regions, corre- 
sponding to different modes n — (dashed gray line), n — 2 (solid black line) as well as instability 
thresholds of traveling wave (TW) bifurcation (dotted blue curves) and Hopf (HH) bifurcation 
(dashed-dotted red curves). The vertical line at a = 1 indicated the onset of spontaneous motion 
(SM). 
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Figure 3. Numerical solutions of Eq. (pQ), obtained for fixed a = 1.05 and different values of the 
delay time r. Numerical simulations have been performed on the two-dimensional square domain 
100 x 100 using a pseudo-spectral method with 512 x 512 grid points, whereas a Runge-Kutta 4 
scheme is employed for the time stepping. Parameters are: a\ — —2, (i2 = —4/3, C = 1, Y = —0.4. 
White arrows indicate the time evolution direction. 



solution of Eq. (|4]), calculated for different a and represents the instability thresholds for 
critical modes as well as thresholds for HH and TW bifurcations. One can see that for a = 1 
both solutions q and qh become unstable. However, for a > 1 the instability type depends 
on the value of r. This is illustrated in Fig. [3j where results of numerical simulations of 
Eq. (pQ) for a = 1.05 and different values of r are shown. For small r all modes are excited 
(see Fig. [3] for r = 0.04) and a stationary bright LS transforms into a breathing dark LS. 
Here, the steady state transfers to a new oscillatory state of a high positive amplitude, i.e., 
a transition from the bright LS to a dark oscillon is observed. If r is increased, qh becomes 



stable, but the spatial modes remain unstable and a soliton ring (r = 0.16) or a labyrinth 
(r = 0.24) is formed. Indeed, a LS moves with a rather large velocity, its shape changes 
(n = 0) and deforms (n — 2). Finally, the impact of the curvature instability results in 
the formation of the labyrinth pattern. Note that labyrinths are not solutions of the SHE. 
However, as mentioned above, there exist parameter regions, where the DSHE behaves like 
a reaction-diffusion system, where the formation of labyrinths is well understood [j^j]. This 
indicates that for certain values of r the DSHE to some extend operates in the regime 
of reaction-diffusion systems. With increasing of r the amplitudes of the unstable modes 
decrease and instead of labyrinths, a moving wave segment emerges on the same time scale 
(r = 0.32). For r = 0.38, both modes n = and n = 2 are stable, i.e., solutions in 
the form of moving LSs are expected. However, close to the instability threshold the real 
parts of corresponding eigenvalues are rather small, what affects the LS's shape on a rather 
large time scale. Finally, a motion with a constant velocity is observed as shown in Fig. [3] 
for r = 0.5. Here, the dynamics of the LS is little affected by critical modes and the LS 
propagates without pronounced change of its shape. 

In order to investigate the stability of the LS in the presence of the TDF from the point 
of view of bifurcation theory we perform the following ansatz 

q = e- R W- v ^q (x)+w(x,t)). (5) 

Here, the operator e _R W' v generates a shift of the pattern to the position R(i) and w(x, t) = 
^j(x) e Xjt is the shape deformation, given by a sum of all stable modes of the delayed 
problem. Our goal is to decompose Eq. © in such a way, that evolution equations for the 
shift R(t) as well as for the change of the shape in terms of w(x, t) are derived. An equation 
R(i) can be obtained by requiring that the shape deformation is orthogonal to the Goldstone 
mode cp^. This condition leads to 

- R(l - (<p£\w)) = -< <^|£'(q ) w) - ( ^|N(w)> (6) 
-a{ v>x|M(x, t, R t - R*_ T , w t , w t _ r )) . 

Here, we introduce the nomenclature Rj := R(t), R- T := R(£ — t), w t := w(x, t) and 
Wj_ T := w(x, t — t), the operator 

N ( w ) = \\ £"(qo)ww + 1 £'"(q )www 
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contains nonlinear contributions, whereas the delay terms are lumped into the operator 

M = E (1 - e ( R *- R ^> v ) (q ( x ) + w(x, t)) 
+Ee (Rt-Rt-.)-v ( Wt _ Wt _ r ) . 

The evolution equation for w(x, t) is obtained by inserting Eq. ([6]) into Eq. ([2]): 

S\V — . ()w — — — 

— = £'w + R— + N + aM . (7) 
ot ax 

Here, the operator £' = SJ — ^ ( ^xl^)/((^x) 2 ) an d other tilded operators are defined 
according to the same rule. Equations ©-(ED show that the deformation of the LS is given 
by a delayed equation and can be treated by a suitable bifurcation analysis. To this end, 
the form of the shape deformation should be specified. As a first approximation, we can 
perform a mode truncation, representing w(x,t) as a superposition of the modes of the 
discrete spectrum 

w(x, t) = u (t)<po(x) + u 2 (*)<p 2 (x)e i2 * + c.c. , 

where the slow functions of time u (t) and u 2 (t) are the amplitudes of modes <£>o( x ) and 
<^ 2 (x). Projection of Eq. ([7]) onto these modes leads to a set of delayed equations: 

uq = fi UQ + a ul + 6 |^2| 2 + c ^o + do u \u 2 \ 2 

+a(u 0it - u 0i t- T ) , (8) 

7i 2 = l^ 2 U 2 + &2 ^2 + C 2 ulu 2 + d 2 \u 2 \ 2 U 2 
+a(u 2 , t ~ U 2 , t -r) • 

Here /x and \i 2 are the eigenvalues of £'(qo), corresponding to the modes n = 0, n = 2, 
and all coefficients are determined by the integrals a = (<^o|£'Vo)/2, &o = (^ol^Vi)? 
c = <<A)|#'Vo>/6, do = (Vo|£"Vo<^>, &2 = (rf'Vo^), c 2 = ((^ |£ / Vo^2)/2, rf 2 = 
(y) 2 |£ /// (^2)/2. Moreover, the delayed equation for R(t) reads 

R = a M - ^ £>*(u M - U kit -r)j (R* " Rt-r) 

-/3|R,-R,_ r | 2 (R,-R,_ r ), (9) 



with ft = ((V^) 2 )/((^) 2 )- Thereby, the elements of the matrix D k are defined as D\- = 
(tpf \dujc(t)/dxj) . A set of nonlinear delayed differential equations (IE])-© is a system of 
order parameter equations, describing the dynamical behavior of a single LS in the vicinity 
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Figure 4. (a) Time evolution of the amplitudes uo(t) (gray) and U2(t) (black), as solutions of 
Eqs. dHJ), calculated for r = 0.38 and a = 1.05. Other parameters are the same as in Fig. [3l (b) 
Dependence of the instability increment on r for three different values of a as result of reduce 
model (|8j)-(j9j) (filled blue triangles) and numerical simulations of Eq. (pQ) (open red squares). 



of the bifurcation point a = 1. Numerical solution of (Ej)-([9j) gives a time evolution of 
the amplitudes of the modes u {t) and U2(t) as well as the evolution of the position R(i). 
In Fig. H] (a) numerical solution of Eq. (|8j), showing relaxation oscillations of both critical 
modes is presented (compare with Fig. [3] for r = 0.38). Note that Eqs. (ED admit solutions 
with u 2 = 0, represented by interplay between the stable decaying breathing mode n = 
and spontaneous motion. Solving the system (|8j)-(|9]) with u 2 = numerically, an instability 
increment, calculated from the time evolution of the amplitude u (£), can be calculated and 
compared with a decay rate, obtained from numerical simulations of Eq. (pQ) (see Fig. H] (b)). 
Now consider solutions with u = 0, corresponding to the case of a spontaneous motion. Let 
us define the drift velocity V(t) := R(t). Close to the bifurcation point a = 1 the velocity is 
slowly varying quantity. Hence, R(t - r) = R(t) - V(t) r + r 2 V(i)/2 + O (d 2 \/dt 2 ) . The 
result is the normal form 



R = V, 

^V = (a-l)V-^/3V 2 V, 



(10) 



which is a normal form of the drift-bifurcation js]. However, since Uq^ = 0, the delay-induced 
drift to the lowest order occurs without change of shape (compare to Fig.[3]for r = 0.5). The 
stationary drift velocity can be calculated directly as V 
order is given by the same relation as calculated in [8|] . 



r 



6 ^ ^ and to the lowest 

p a 



To conclude we have shown that a SHE subjected to TDF exhibits complex spatio- 
temporal dynamics. The existence of labyrinths indicates that in certain delays the DSHE 
to some extend behaves like a reaction-diffusion system. The mode truncation, neglecting 



contributions of the strongly damped modes of the continuous spectrum leads to a system of 
order parameter equations, describing the behavior of the single LS in the DSHE. Moreover, a 
normal form of the delay-induced drift-bifurcation is obtained, showing that induced motion 
to the lowest order arises without change of the shape. In addition, presented analytical 
results are derived in general form and can be extended to other spatial extended systems. 

We thank Andrei G. Vladimirov for fruitful discussions on the topic. 
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